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ABSTRACT 

Cosmic rays with energies exceeding 10^° eV have been detected. The origin of these 
highest energy cosmic rays remains unknown. Estabhshed astrophysical acceleration mech- 
anisms encounter severe difficulties in accelerating particles to these energies. Alternative 
scenarios where these particles are created by the decay of cosmic topological defects have 
been suggested in literature. In this paper we study the possibility of producing the high- 
est energy cosmic rays through a process that involves formation of metastable magnetic 
monopole-antimonopole bound states and their subsequent collapse. The annihilation of the 
heavy monopole-antimonopole pairs constituting the monopolonia can produce energetic nu- 
cleons, gamma rays and neutrinos whose expected flux we estimate and discuss in relation 
to experimental data so far available. The monopoles we consider are the ones that could 
be produced in the early universe during a phase transition at the grand unification energy 
scale. We find that observable cosmic ray fiuxes can be produced with monopole abundances 
compatible with present bounds. 
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1 Introduction 



The physics and astrophysics of cosmic rays (CR) of ultrahigh energy (UHE) (i.e. with energy 
above about 10^® eV) constitute a subject of much intense research [1, 2, 3] in recent times 
both in terms of new experiments as well as new theories. UHE CR with energies exceeding 
10^° eV have been detected [4, 5, 6, 7, 8, 9, 10]. The Haverah Park experiment [4] reported 
several events with energies near and slightly above 10^° eV. The world's highest energy CR 
event detected recently by the Fly's Eye experiment [7, 8] has an energy ~ 3 X 10^° eV. The 
event of energy ~ 1.1 X 10^° eV recorded by the Yakutsk experiment [5, 6] has almost the same 
arrival direction as that of the Ely's Eye event. More recently, the AG AS A experiment [9] 
has also reported an event with energy (1.7-2.6) X 10^° eV [10]. 

The existence of UHE CR, especially, the highest energy cosmic rays (HECR) (i.e. with 
energy above 10^° eV) poses serious challenge [11, 12, 13, 14, 15, 16] for conventional astro- 
physical acceleration mechanisms [17] that attempt to explain the origin of these particles 
in terms of acceleration in special astrophysical sites like supernova shocks, pulsar mag- 
netospheres, galactic wind termination shocks or relativistic shocks associated with active 
galactic nuclei and radio galaxies. In this last case acceleration up to around 10^^ eV seems 
to be possible by stretching the "reasonable" values for the shock size and the magnetic field 
strength at the shock somewhat [18]. However, at least for the highest energy Ely's Eye and 
Yakutsk events mentioned above, as also for the more recently detected AGASA event [10], 
there seem to be no suitable extragalactic objects such as AGNs or rich galaxy clusters near 
the observed arrival directions and within a maximum distance of about 50Mpc, this upper 
limit on the possible source distance being set by considerations of energy loss during prop- 
agation [13, 14, 15, 16]. Thus it is difficult to associate these highest energy events with any 
known astrophysical sources. 

The difficulties encountered by conventional acceleration mechanisms in accelerating par- 
ticles to the highest observed energies have motivated recent suggestions [19, 20, 21, 15] that 
the underlying production mechanism of HECR could instead be of a non-acceleration na- 
ture, namely the decay of supermassive elementary "X" particles related to Grand Unified 
Theories (GUTs). Sources of such particles today could be topological defects (TDs) [22] 
formed in the early universe during phase transitions associated with spontaneous breaking 
of symmetries implemented in these GUTs. This is because TDs like cosmic strings, domain 
walls, superconducting cosmic strings and magnetic monopoles are topologically stable but 
nevertheless can release part of their energy in the form of these X particles due to physical 
processes like collapse or annihilation. The X particles, with masses of the typical GUT 
scale which is generally much higher than 10^^ eV, decay into leptons and quarks, the latter 
ones finally hadronizing into jets of hadrons and giving rise to HECRs. In this scenario, the 
observed HECR are due to collapse or annihilation of TDs at relatively "close" distances 
(< 50Mpc) from earth. 

The predicted spectra of UHE particles in the TD scenario are determined essentially 
by the physics of hadronization of quarks, i.e., by Quantum Chromodynamics (QCD). In 
this sense the shapes of the spectra are universal (i.e., independent of the specific process 
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involving any specific kind of TDs), especially at the highest energies where cosmological 
evolutionary effects are negligible (except for neutrinos). The overall contribution of different 
TD processes to the UHE cosmic ray flux would, however, be different. 

Among the various kinds of possible TDs, the case of cosmic strings is perhaps the one 
that has been studied most extensively in terms of their formation and subsequent evolution, 
principally because they provide an attractive theory of formation of galaxies and large-scale 
structure in the universe. It was, therefore, natural to investigate first the various possible 
UHE CR producing processes involving cosmic strings. Almost all the mechanisms studied 
so far in this connection [20, 23, 24] involve closed loops of cosmic strings. It turns out that 
cosmic strings can give measurable contribution to UHE CR fiux only if there is a mechanism 
by which a fraction ~ I0~^ of the energy of all closed loops of cosmic strings chopped off the 
long (i.e., larger than horizon size) segments of strings at any time is dissipated in the form 
of X particles on a time scale much smaller than the time scale of energy loss of these loops 
through gravitational radiation. One process which would satisfy this criterion involves the 
"completely" collapsing class of loops [20] or loops which undergo multiple self-intersections 
and break up into a large number of small loops (rather than a small number of large loops) 
within one period of oscillation of these loops. These loops are, however, rather special in 
the sense that they have to be fine-tuned to one of the collapsing configurations to within 
a length scale of the order of the width of the string, which is a microscopic scale of order 
~ I0~^^ cm for GUT scale cosmic strings. Gill and Kibble [25] have recently argued that 
these "smooth" loops are unlikely to be formed in any significant number except possibly 
at very early times. To explain HECR, however, the collapse of the loops must occur in 
relatively recent epochs. Gill and Kibble have, therefore, argued that processes involving 
ordinary cosmic strings are unlikely to yield a measurable fiux of UHE CR. It is, however, 
conceivable [26] that hybrid systems of TDs such as light domain walls bounded by GUT 
scale cosmic string loops could form in the early universe; the domain walls in these systems 
could aid the complete collapse of the string loops that form the boundaries of these walls. 
It is also conceivable [27] that gravitational radiation back reaction effects may smooth out 
higher frequency wiggles leaving only the lowest-frequency mode waves on the string loops 
making them collapse completely [28] such that a significant amount of UHE CR could be 
produced. These possibilities, however, remain to be studied in detail. 

Another potential TD-source of HECR are the saturated superconducting cosmic string 
(SCS) loops [29, 21]. If and when an SCS loop achieves a certain saturation current [29] 
the massive charge carriers that carry the electric current on the string are ejected from 
the string. The subsequent decay of these massive charge carriers can produce UHE parti- 
cles [21]. There are, however, several issues in this context that remain to be settled. The 
cosmological evolution of current-carrying SCS is much more complicated than evolution of 
"ordinary" cosmic strings. In particular, it is possible [30] that SCS loops may never achieve 
the saturation current at all. Another issue of debate [31, 21] is whether or not the charge 
carriers, assuming they are ejected from the string, can get out of the immediate vicinity 
of the string before decaying. The decay of the massive charge carriers must occur outside 
the region of strong magnetic fields surrounding the strings; otherwise, the energetic decay 



2 



products would rapidly lose energy by interacting with the strong magnetic field and so 
would not survive as UHE particles. 

These issues certainly remain to be studied in detail. At the same time, it is worthwhile, 
in our opinion, to study the possibility of production of measurable fiuxes of UHE cosmic 
rays by other possible processes involving other kinds of TDs. With this motivation, we 
investigate in this paper another possible TD-process of generating UHE cosmic rays, namely 
the collapse of relic magnetic monopole-antimonopole bound states (monopolonium), first 
suggested by Hill [32]. 

Magnetic monopoles (simply "monopoles" hereafter) are one kind of TD solutions in 
spontaneously broken non-abelian gauge theories that are allowed in essentially all GUT 
models. If nix denotes the mass of a typical gauge boson in a GUT, then the monopole 
mass ruM is given by tum ~ ctx^mx, where ax is the dimensionless "unified" gauge cou- 
pling strength at the energy scale nix in the GUT model under consideration. The core 
of a "GUT monopole" has a radius ~ m^^. Eormation of massive monopoles in the early 
universe, through Kibble mechanism [22], is inevitable in most GUT models, and leads to 
the well-known "monopole problem" (see, for a review, Ref. [33]). The "problem", to re- 
mind the reader, refers to the "standard model" prediction of relic abundance of monopoles 
which, on various empirical grounds, is unacceptably large. (Here, by "standard model" 
prediction we mean prediction made within the context of the simplest GUT models and 
the standard Big-Bang cosmology). Several mechanisms including cosmic infiation [34] have 
been proposed for reducing the relic monopole abundance to acceptable levels (see, for a 
review, [33]). Eor instance, interesting relic abundances could have been produced thermally 
during reheating after the universe has gone through an infiationary phase. However, we will 
not discuss those mechanisms in this paper. Instead, we will simply assume that monopoles 
exist in the universe at a level of abundance compatible with known experimental [35] and 
phenomenological [36, 37] upper bounds. 

Since monopoles are topologically stable, the only way of getting rid of monopoles (Ms) is 
to make them annihilate with antimonopoles (Ms). The "standard" way (see, e.g., Refs. [38, 
39]) of achieving this relies upon mechanisms to capture M — M pairs in metastable bound 
states which spiral in and finally collapse resulting in annihilation of the Ms and the Ms 
that have been captured in bound states. However, as is well known, this typically is a 
slow process and fails in "solving" the monopole problem, i.e., the Ms and Ms do not 
annihilate fast and early enough for the Universe to avoid being monopole dominated. //, 
however, the monopole problem is "solved" to start with (by some mechanism which we do 
not concern ourselves with in this paper), i.e., if the universe is never monopole dominated 
to start with, which we shall assume to be the case in this paper, then the late annihilation 
of the monopoles is, in fact, precisely the mechanism that we need from the point of view of 
generating UHE cosmic rays which, as we know, must be produced only in the contemporary 
cosmic epoch (i.e., at low redshifts); the potential UHE cosmic ray particles resulting from 
M — M annihilation occuring in the earlier epochs would thermalize by interacting with the 
dense background medium and hence would not survive as UHE particles. So monopoles 
formed at a GUT-scale phase transition in the early universe and annihilating in the recent 
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epochs provide us with an attractive scenario of production of UHE cosmic rays provided 
they exist in the Universe in sufficient numbers. The main aim of the present paper is, in 
fact, to try to estimate the monopole abundance required in order to generate enough UHE 
CR flux as observed. Furthermore, we discuss the predicted particle spectrum in relation to 
experimental data so far available. 

We organize the rest of the paper as follows: In Sec. 2 we describe the relevant frame- 
work of monopolonium physics. In Sec. 3 we calculate the contribution of monopolonium 
collapse to the UHE CR flux in terms of the model parameters and estimate the abun- 
dance of monopolonia (relative to that of monopoles at the relevant time of formation of 
the monopolonia) required to produce enough flux of the highest energy cosmic rays as ob- 
served. In section 4, we then study, following the original analysis of Hill [32], whether the 
required monopolonia abundance established in section 3 are realizable. Einally, in Sec. 5 
we summarize our findings. 

Wherever appropriate we shall use natural units with h = c = ks = 1^ ks being the 
Boltzmann constant. Quite often, however, we keep one or more of these quantities explicitly 
in some formulae for convenience, which should be clear from the context. Also, we assume 
a "fiat" (Oo = 1) universe with a Hubble constant Hq = lOO/ikmsec"^ Mpc~^ with h=0.75 
taken in the actual numerical calculations. 



2 Monopolonium Physics: A Brief Introduction 

Monopolonium [32] is a possible bound state of a magnetic monopole and an antimonopole. 
Let us first consider, classically, a monopole (Af) and an antimonopole (Af) separated by a 
distance r and bound in a circular orbit around each other. The non-relativistic energy of 
the system can be written as 

E = -niMto'^r'^ - gl,lr , (1) 

and since 

— = r , (2) 

r 

we have 

E = -\gljT, (3) 

where rriM = rnml'^ is the reduced mass, u is the orbital angular frequency, and is the 
magnetic charge of the monopole. The Dirac quantization condition relating the magnetic 
charge g^ and electric charge e is (keeping fi and c) [40, 41] 

egm_^ . . 

he 2 ' ^ ' 

where the integer N is the "monopole number". Eor = 1 monopoles, the "magnetic" 
fine-structure constant is given by 

gl I I m 
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where a^. = j ific) ~ 1/137 is the "electric" fine-structure constant. We will ignore the 
"running" of the coupling strengths and 

"Classical" monopolonium is, of course, unstable (like the "classical" hydrogen atom). 
A Bohr model of monopolonium can, however, be constructed. In the Bohr model, not all 
values of r are allowed. Monopolonium can exist only in certain discrete states characterized 
by the principal quantum number n given by 



n^a^ , (6) 



where n = 1, 2, 3, . . ., and 



= = 8ae = Sttettx (7) 

mmgin \mMcJ \mxc J 

is the Bohr radius of the monopolonium state. In writing Eq. (7) we have used Eq. (5). Now 
Eq. (3) can be written as 

E = -^. (8) 
where is the effective "Rydberg constant" for monopolonium, and is given by 

= IrhM^f ^ 2.5 X 10^^ al {j^^) GeV . (9) 

(Compare with this the Rydberg constant for the hydrogen atom, 13.6 eV). 

Note from Eq. (7) that ^ fi/{mMc) < fi/{mxc). In other words, the monopolonium 
ground state (n = 1) is one in which the cores of the monopole and the antimonopole overlap 
strongly. Clearly, then, monopolonium does not really exist in the n = 1 state because the 
M and the M would annihilate each other. Presumably, when monopolonium is formed 
in the early universe, the M and the M capture each other in a state with n ^ 1. This 
state then undergoes a series of transitions to tighter and tighter bound states by emitting 
initially photons and subsequently gluons, Z bosons, and finally the GUT X bosons. At some 
stage during the collapsing process, the cores of the M and the M begin to overlap (when 

1/2 

r ~ 2/mx, i.e., when n reduces to ~ [rriM / {4:C(emx)] ) and subsequently annihilate each 
other. When this happens, the bound state is destroyed and the energy contained in the 
system is released in the form of various particles. We are interested in estimating the fiux 
of UHE CR produced by these final M — M annihilation events associated with collapsing 
monopolonia in the universe. Actually, the gluons emitted by Larmor radiation during the 
collapsing process as well as the quarks from the decay of Z-bosons radiated by monopolonia 
will also hadronize and produce energetic hadrons which will contribute to the total CR 
fiux. We will, however, restrict ourselves here to estimating the CR fiux due to the final 
M — M annihilation events discussed above. Therefore, what we will estimate here will be a 
lower limit to the total possible contribution to the UHE CR fiux from monopolonia in the 
universe. 
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The monopolonium lifetime r can be calculated [32] by using the dipole radiation formula 



where Eitj) = — (/^/(2rj), rj being the radius of the monopolonium at the time of its 
formation, and E{tc) = —g'^/{2rx) with rx ~ itt-x^- 

In the situation relevant for our case, rj ^ rx ^ 10~^^cm. The lifetime r = tc — tj is, 
therefore, given by 



Note that r oc rj. Thus, for example, r ~ 40 days if rj ~ Ifm, whereas r ~ 10^^ yr 
if rj ~ 1 nm. In other words, depending on the initial radius of the bound state at its 
formation, some of the monopolonia formed in the early universe could be surviving in the 
universe today and some would have collapsed in the recent epochs. 

An operational definition of "formation" of monopolonium states in the universe can 
be taken as follows. At any time t when the temperature of the universe is T, Af — Af 
bound states with binding energy Ej, = rjT (where rj > 0(1) is an unknown parameter at 
this stage) "freeze out" or "form" and start to collapse. The parameter rj incorporates the 
requirement that the formed monopolonia not be thermally dissociated. (We are, of course, 
implicitly assuming here that monopoles interact efficiently enough with the background 
plasma of thermal electrons and photons at the relevant times of formation of the bound 
states; see below). The abundance of monopolonium formed can in principle be estimated 
if a specific mechanism of formation is given. We will not go into the discussion of any 
specific monopolonium formation mechanism here. Instead, we will first try to estimate, on 
a phenomenological basis, the monopolonium abundance (relative to a given abundance of 
monopoles) required for the scenario to yield measurable fiux of UHE cosmic rays. In section 
4, we will estimate the relative monopolonium abundance within the framework of a Saha 
formalism applied to a system of Ms, Ms and M — Ms in thermal equilibrium, and compare 
the abundance so obtained with the phenomenologically required estimate. We shall see 
that within the framework of the thermal equilibrium analysis the relative monopolonium 
abundance is essentially parametrized by the parameter rj mentioned above. 

We are interested in monopolonium states that are collapsing in the present epoch, i.e., 
tc ^ ^eg, where teq is the time of equality of radiation and matter energy densities. On the 
other hand, with the operational definition of formation of monopolonia mentioned above 
and the life-time formula, Eq. (12), we have, tj <^ tpq, and, therefore. 




T = 






(13) 
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where 

Tj = 1.56 X 10-3 C'^' ( —) GeV . (14) 

V sec/ 

In Eq. (14) is the total number of effective relativistic degrees of freedom determining 
the energy density of radiation in a radiation dominated universe [33]. Eq. (13) reveals 
that monopolonium states undergoing final collapse today must have been formed with 
a binding energy E^j ~ 0.5 (mM/10^^ GeV)^''"^ MeV, corresponding to a temperature of 
the universe at the time of formation, Tj, also of this order, i.e., around the epoch of 
primordial nucleosynthesis. Comparing the typical time scale for monopole-plasma energy 
exchange [38], r, ~ 6.58 x lO'^ 5*"^^^ (mw/lO^^ GeV) ( MeV/T)^ sec, with the Hubble time 
(age), tu ~ 2.42 (/* ^''^ ( MeV/r)^ sec, of the radiation (i.e., relativistic plasma) dominated 
universe, we see that < t„ as long as niM < 3.7 X 10^^ GeV. In other words, the monopoles 
may be assumed to be in thermal equilibrium with the electron-positron-photon plasma at the 
time of formation of the monopolonium states relevant for our considerations. However, the 
e'^e~ annihilations at T ~ 0.3 MeV [33] significantly reduces the effectiveness of monopole- 
plasma scatterings in maintaining thermal equilibrium of the monopoles. Thus, while the 
relevant bound states may be assumed to be formed when the monopoles are still in thermal 
equilibrium (although only marginally so if we take into account the requirement that rj > 
1 so that Tf < -Eft/), their subsequent "spiraling in" and collapse essentially occurs in a 
situation in which the monopoles are effectively decoupled from the background medium. 
This justifies, albeit a posteriori^ our use of the "vacuum" dipole radiation formula, Eq. (10), 
in calculating the lifetime of the relevant monopolonia at least at the level of approximation 
adopted in this paper. 

We shall use Eqs. (13) and (14) for calculating tj for a given value of tc or vice versa. 



3 Contribution of Monopolonium Collapse to UHE 
CR 

3.1 The Rate of X Particle Production 

The number of X particles released due to annihilation of M and M constituting monopolo- 
nium can be simply taken to be 2mM I mx ■ The number density of X particles released per 
unit time, dnx /dti^ due to collapsing monopolonia in the universe at any time ti = tc is then 
given by 

rln V rln': '?'m 

(15) 



dnx _ dnl^j^ 2m m 



dti dti nT-x 

where 



dti dEb f dti \^ ~^ Zf , 



(16) 



Here f^^jg is the number density of monopolonia collapsing at the time ti and f^^^^g is their 
density at formation. Also Zi and zj are the redshifts corresponding to the times ti and 
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tf. The redshift factor in Eq. (16) takes care of the dilution of the monopolonium number 
density due to expansion of the universe. 

3.2 Injection Spectra of Nucleons, Gamma Rays and Neutrinos 

We now assume that each X-particle decays into a lepton and a quark each of approximate 
energy mx/2. The quark hadronizes by jet fragmentation and produces nucleons, gamma 
rays and neutrinos, the latter two from the decay of neutral and charged pions in the hadronic 
jets. The lepton can also generate further particles by interacting with the background 
medium. But we expect the hadronic route will generate by far the largest number of 
particles, and we will concentrate on these. 

The spectra of the hadrons in a jet produced by the quark are, in principle, given by 
QCD. Suitably parametrized QCD motivated hadronic spectra that fit well the data in 
collider experiments in the GeV-TeV energies have been suggested in the literature [32]. 
Below, we shall illustrate our results by using the QCD motivated spectra suggested by 
Hill [32]. It is to be kept in mind, however, that there is a great deal of uncertainty involved 
in extrapolating the formula that describe the "low" energy data to the extremely high 
energies as in the present situation. To study the sensitivity of our results to the assumed 
hadronization spectrum, we will also present the results for injection spectra suggested in 
Ref. [42] on certain phenomenological grounds. 

The injection spectrum, i.e., the number density of particles produced per unit time per 
unit energy interval, for the species a= nucleons (N), gamma rays (7) and neutrinos (;/), can 
be written as 



where x = 2Ei/mx^ Ei being the energy at injection and dNa/dx is the effective fragmenta- 
tion function describing the production of the particles of species a from the original quark. 
We will consider the following two cases for the fragmentation function: 
(1) QCD-motivated injection spectra 

In this case, the total hadronic fragmentation spectrum dN^/dx is taken in the form [32] 



where the lower cutoff Xq is typically taken to correspond to a cut-off energy ~ 1 GeV. 
Assuming a nucleon content of ~ 3% and the rest pions, we can write the fragmentation 
spectra as [19, 43] 




(17) 




^ if Xq < X < 1 



otherwise 





(19) 



dN^[x) 



( 




1 1 dNhix') 



dx 



X X tlx 
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The neutrino spectrum above includes only the (z/^ + t7t)'s resulting from the first-stage 
of the charged pion decay, i.e., from tt^ — )• //^ + (z?^). The further decay of the /i^, i.e., 
/i^ — )• + Ueiue) + i^iiii^ij.), produces additional neutrinos of approximately the same spectra 
as those of the "first stage" neutrinos. Thus altogether we may expect to have roughly twice 
as many muon neutrinos as given by the last line of Eq. (19). For a conservative estimate, 
however, we shall show our results only for the "first stage" neutrinos mentioned above. 

(2) Phenomenological injection spectra 

Recently, injection spectra somewhat different from the QCD motivated injection spectra 
described above have been suggested by Chi et al [42] on the following phenomenological 
grounds. UHE gamma rays as well as protons generate lower energy gamma rays by 7 — 75 
and p — 76 collisions with the photons (75) of the background radiation fields. The elec- 
tromagnetic component of the energy lost by the photons and protons in these collisions 
cascades down to lower energies by electromagnetic cascading in the universal radio back- 
ground (URB), the cosmic microwave background (CMBR) and in the infrared background 
(IRB) (in order of decreasing energy of the propagating photon). The measured fiux of ex- 
tragalactic gamma rays in the 100 MeV energy region [44] provides constraints on the form 
of the nucleon and gamma ray injection spectra at energies above ~ 5 X 10^^ eV. Based on 
these considerations, Chi et al [42] have suggested injection spectra which we describe by 
the following fragmentation functions: 



dNN[x) 



dx 



Anx-^-% (20) 



dK 



X 



dx 
dNjx] 



X 



A}y X 



with 



anc 



dx 

A„ = A.,^ 0.01 , (211 



A^ /lO^^GeV'"-' 



, _ 0.028 . 22 

An \ mx J 

The condition (22) comes from the requirement [42] that the ratio of photon-to-nucleon 
{^f/N) at injection at energy E = 10^° eV, i.e., at x = 2 X 10^°eV/mx be ~ 60. The spectra 
(20) are assumed to be valid above a suitable lower cutoff Xq < 10""*. 



3.3 The Evolved Spectra 

The evolution of the spectra is governed by energy-loss and/or absorption of the particles 
as they propagate through the extragalactic medium. We are interested in calculating the 
expected diffuse fiux assuming that the monopolonia are distributed uniformly. The general 
expression for the expected diffuse fiux today [t = to), i.e., the number of particles crossing 
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per unit area per unit time per unit solid angle per unit energy interval at an energy Eq^ can 
be written as [20, 23] 



j{Eo) = ^cto dz,{l + ^0 '^ $(^., z,) , (23) 

OTT Jo uEq 

where Zi is the injection redshift corresponding to the injection time t^, Ei[Eo^Zi) is the 
necessary injection energy, and the maximum injection redshift Zi^max{Eo) is determined 
from the condition Ei{zi^raax^ Eq) < mx/2. 

We will use the continuous energy-loss approximation [45, 46] for all particles. The 
general energy-loss equation in terms of redshift z can be written as [20, 23, 45] 

^ 1 + z i^o 

The first term is due to redshift of the (relativistic) particle energy and the second term 
describes losses due to interactions with the background medium in terms of the inverse 
energy-loss time scale /3{E) = —[dE/dt)/E. For a particle observed at the earth with an 
energy Eq the necessary injection energy Ei[zi) at redshift Zi can be found by integrating 
Eq. (24) back from z = to z = Zi with Ei[z = 0) = Eq. 

For nucleons we shall use the expression for /3{E) derived in Ref. [45]. The important 
process here is the Greisen-Zatsepin-Kuzrhin (GZK) [47] effect in which UHF nucleons above 
a threshold energy (~ 6 X 10^^ eV) lose energy drastically due to photopion production off 
the CMBR photons. This gives rise to the onset of a sharp fall ("cutoff") of the extragalactic 
UHF nucleon fiux. The interactions of UHF nucleons with the CMBR also cause secondary 
gamma rays and neutrinos. We shall neglect their contribution to the total gamma ray and 
neutrino fiux calculated below thus getting lower limits for the fiuxes of these particles. 

The evolution of the UHF gamma ray spectrum is mainly governed by absorption of the 
UHF photons through e"'"e~ pair production on the CMBR and on the URB. Actually, under 
certain circumstances, the propagating photons give rise to electromagnetic cascades [43]. 
This results in an increase of the effective penetration length of the UHF gamma rays, 
which, in turn, has the effect of increasing the final gamma ray fiux. The cascading effect, 
however, depends rather strongly on the strength of the intergalactic magnetic field which is 
rather uncertain. We will ignore here the cascading effect and consider only the absorption 
of the UHF photons on the CMBR and URB photons. This will give us a conservative 
estimate of the final gamma ray fiux leading to a conservative estimate of the required relative 
monopolonium abundance, i.e., the actual required monopolonium abundance should be even 
lower than what we estimate. We will take the absorption lengths for UHF gamma rays as 
given in Ref. [43]. The relevant absorption lengths being small compared to the Hubble size 
of the universe, the cosmological redshift term in Fq. (24) is essentially immaterial. Thus 
for gamma rays we consider only absorption and no energy-loss of the propagating gamma 
rays. 

For UHF neutrinos the dominant process relevant for evolution of the spectrum is the 
absorption [48, 49] of neutrinos through the process p -\- Pj, ^ ff^ where / = e, r, u, J, s, c. 
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etc., and Pj, represents the thermal background neutrinos (which have a present temperature 
~ 1.9 K). It has been pointed out [50] that UHE neutrinos would also generate a neutrino 
cascade effectively increasing the "neutrino horizon" of the universe, which, in turn, has 
the effect of increasing the final fiux of neutrinos. For example, for an X-particle mass 
nix = 10^^ GeV the fiux at 10^^ eV is increased by about a factor 6 [50]. Again, since we 
are interested in a conservative estimate of the fiux, we will ignore this cascading effect here. 
We will take the absorption redshift for UHE neutrinos as given in Ref. [19]. 

3.4 Required Monopolonium Abundance: An Approximate An- 
alytical Estimate 

With the knowledge of energy attenuation lengths and absorption lengths for various par- 
ticles as described above, we can calculate the expected fiux of particles at earth. We will 
discuss the results in detail in the next subsection. Here we present an approximate ana- 
lytical estimate of the required monopolonium abundance (relative to a given abundance of 
monopoles at the relevant time of monopolonium formation) that would yield sufficient fiux 
of UHE cosmic rays as observed. We will do this by estimating the expected UHE gamma 
ray fiux and matching it with the observed fiux at a given energy. 

We define the relative monopolonium abundance at formation, i^j, for the monopolonia 
collapsing in the present epoch, as 



where the integral is understood to be taken over a range of binding energies at formation 
for which the corresponding monopolonia would be collapsing within, say, one Hubble time 
period (~ tc) at the time of collapse tc which for the CR fiux today is to, the present age of 
the universe. 

Neglecting cosmological effects, the expected UHE gamma ray fiux, j(£'o), at an observed 
energy Eq is simply 



where A(£'o) is the absorption path length. Eor illustration, we shall use the QCD motivated 
gamma ray injection spectrum described by Eqs. (17), (18) and (19). At an energy Eq for 
which X = 2Eo/mx ^ 1, we get 



Using Eqs. (15) and (16) for dnx /dtc^ and taking the gamma ray path length from Ref. [43], 
we get, for = 3 X lO^OeV, 




(25) 




(26) 




(27) 
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E'^jiEo) ~ 7.3 X 10^^ • m^sec^ 



A \ ^mM\ ( nix ,,2 -2 -1 -1 /oo^ 
I I — — I eV m sec sr . z8 



VlOMpc ; \mxJ VlQis GeV 

Comparing this with the observed flux (corresponding to the Fly's Eye's observed highest 
energy event), which is Elj{EQ = 3 X 10^° eV) ~ 2 X 10^^ eV^ sec~^ sr~^, we get 



^^-2.74x10"^° r^]\—^]\ m-^sec-^ 29 

ate \ X J KtiimJ V ''^x I 



Now, using Eqs. (16), (29) and (25), and noting that 

1 + ^ 

T 

M 



1+ZjJ 



(30) 



where the superscripts c and / refer to the quantities at the time of collapse and formation, 
respectively, we get, taking ti = tc = to = 2.057 X 10^'^h~^ sec, the following expression for 
the required relative monopolonium abundance: 

In Eq. (31) Om is the mass density contributed by monopoles (in units of the critical density of 
the universe), which is related to the ratio of monopole-to-photon number densities, um/it-^^ 
through the relation 



and n^{to) is the photon number density in the universe today. We should perhaps mention 
that throughout this paper we assume the monopolonium formation rate to be low enough 
(which is justifled a posteriori) that the ratio um/it-^ is approximately constant. 

Thus, Eq. (31) gives us a rough estimate of the required monopolonium abundance 
relative to any given monopole abundance. The results of our numerical calculations de- 
scribed in the next subsection conflrm this expectation. As expected we need larger relative 
monopolonium abundance for smaller monopole abundance. On the other hand, for a given 
monopole number density, the larger the value of nix (or equivalently mj^f), the larger is 
the flux produced (see, e.g., Eqs. (27) and (28)), and hence, smaller is the required relative 
monopolonium abundance (note that Om in Eq. (31) is proportional to nix)- 



12 



3.5 Results 



We now describe the results ol our lull numerical calculations of the flux and the required 
relative monopolonium abundance. 

The recent results of the Fly's Eye experiment [7] indicate a change of composition of the 
UHE cosmic rays from a predominantly heavy component at energies below about 2 X 10^® eV 
to a predominantly light component above that energy. Further, this composition change is 
correlated with a dip in the overall spectrum [7]. The Fly's Fye group obtained a good fit 
to the observed spectrum by superposing a steeper galactic component of iron nuclei and a 
fiatter extragalactic component of protons. If we assume that the extragalactic component is 
dominated by monopole annihilation at the high energy end we can determine the necessary 
relative monopolonium abundance mentioned above by normalizing the calculated fiux to this 
light component. In the Fly's Fye as well as in the AG AS A data this high energy end which 
we shall use as our normalization point is located at about Eq = 10^^''' eV. Between this 
energy and the respective highest energies measured there is a gap in both data sets which 
could indicate a break in the spectrum. It will be clear from Figs. 1 and 2 presented below 
that, because of the specific shapes of the spectra, normalization at a lower energy would 
be inconsistent because then the calculated fiux would exceed the observed fiux at a higher 
energy. In this sense, the above normalization procedure is the only consistent normalization 
in our case. Note also that since gamma rays and nucleons can not be distinguished by the 
Fly's Fye at these energies we have to use for normalization the combined nucleon and 
gamma ray fiux given by our model. 

In our numerical calculations we use = 34.25 from Fq. (5) and ax = 1/30. The 
monopole abundance Om is constrained by various bounds [33]. Currently, the most strin- 
gent (and theoretically well-motivated) constraint on Om is that given by a recent modifica- 
tion [37] of the original Parker bound [36], giving Om/^^ < 3.99 x IQ-^ {niM/lO^^ GeVf for 
monopole velocities ~ 10~^c. (Note that current experimental upper limits [35] on Om are 
still about a factor of 5 weaker than the recent Parker bound). For reference we will therefore 
use the improved Parker bound [37] for Om- The monopole mass is fixed by nix and ax- 
We will show the results for two values of mx, namely, nix = 10^^ GeV and 10^^ GeV. 

Figs. 1 and 2 show the nucleon, gamma ray and neutrino spectra obtained by the nor- 
malization procedure mentioned above. Figs. 1 and 2 correspond to the QCD motivated and 
the phenomenological injection spectra described in section 3.2, respectively. As expected, 
the gamma rays and neutrinos far outnumber the protons at energies above ~ 10^° eV. 

We can compare the predicted and the "measured" fiuxes at the highest energies in the fol- 
lowing somewhat more quantitative way: The integral CR fiux above 1.7 X 10^° eV estimated 
from the highest energy Fly's Fye and AGASA events is ~ 0.7 X 10~^^ m~^ sec~^ sr~^ [51]. 
The corresponding gamma ray fiux expected within our model is about a factor 10 (for QCD 
motivated injection, see Figs. 1) and 3 (for phenomenological injection, see Figs. 2) larger. 
In contrast, the expected nucleon fiuxes are smaller by about a factor 3 and 5, respectively. 
Using the additional information about the event energies and comparing with the predicted 
differential fiux (see error bars in Figs. 1 and 2) our model is at least for the phenomenologi- 
cal injection spectrum in good agreement with gamma rays as primaries for the two highest 
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energy events whereas nucleons seem to be disfavored. 

Taking the enhancement factor due to neutrino cascading mentioned in section 3.3 into 
account the predicted differential neutrino flux is between one and two orders of magnitude 
smaller than current limits on the flux of deeply penetrating particles from the Fly's Eye 
experiment [52, 53]. Furthermore, the predicted integral neutrino flux above ~ 10^° eV is 
smaller than current limits from the Frejus detector [54] by factors of ~ 10~^ (QCD motivated 
injection) and ~ 10"'* (phenomenological injection). 

Finally, Figs. 3 and 4 show the contours of the required values of i^, the relative monopolo- 
nium abundance (or equivalently, the value of the parameter ?y, see section 4 below), obtained 
by the normalization procedure described above, in the Om — plane. The shaded regions 
are excluded by the Parker bound on Om- The sharp upturn of the curves for very low values 
of nix in Figs. 3 and 4 is simply a reflection of the fact that for too low values of nix the 
maximum available energy ~ nix /2 becomes lower than the energy of the cosmic rays under 
consideration. 



4 Equilibrium Estimate of Monopolonium Abundan- 
ce and the Resulting Cosmic Ray Flux 

In this section, we estimate the relative monopolonium abundance in the universe by using 
the classical differential version of the Saha equation applied to a system of Ms, Ms and 
M — M bound states in thermal equilibrium with the background relativistic plasma (see 
the discussions at the end of section 2). Following Hill [32] we shall use the Maxwell- 
Boltzmann distribution for the number density of monopoles to estimate the number density 
of monopolonia. As Hill pointed out, the equilibrium analysis is expected to yield a lower 
limit to the true monopolonium abundance, i.e., consideration of any specific mechanism of 
monopolonium formation is likely to increase the monopolonium abundance, which would in 
the end increase the contribution to the total cosmic ray fiux. In other words, the equilibrium 
monopolonium abundance estimated below probably gives a conservative estimate of the 
possible contribution of monopolonia to the cosmic ray fiux. 

Within the framework of an equilibrium analysis, the operational definition of monopolo- 
nium "formation" adopted in section 2 implies that at a time t when the temperature of the 
universe is T, monopolonia with binding energy Ej, = rjT are "formed" with an abundance 
equal to the equilibrium abundance of M — M bound states of binding energy rjT. The 
equilibrium abundance can be calculated by using the Saha formalism. Thus monopolonia 
formed between times tj and tj -\- dtj have binding energies between E^j and E^j -\- dEbj. 
Their number density can be written as [32] 

where E^j = rjTj^ and um (it-m) denotes the number density of monopoles (antimonopoles) 
at the time tj. Using = 2([3)Tf /tt'^ (C(3) = 1.202) for the photon density at temperature 
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Tf and Eq. (32) we can write Eq. (33) in the form 

J 



^.,9TxlO-^<,r-e..n„;.'(^)-'. ,34, 



We thus see that the quantity rj parametrizes the bound state abundance for a given monopole 
abundance Om- The relative monopolonium abundance at formation, i^/, is then given by 
Eq. (25). 

Using Eqs. (15), (16) and (34), and converting redshifts to time by using standard cos- 
mological relations [33], we can write 

dnx 4C(3)2 g 3 fmM\ e"^ ( ^mV ,2,-^ 



dti 37r ™ \mx / (27r?y)'^/^ \ n 



-57 3 (o .2Yf \-\-3l2,n 



W^yA ( To \' fto^' 



to / V2.75K 



-3 „„„-! 



) ^^^^ ' 



where To is the photon temperature today and to is the present age of the universe. 
Introducing the generic form [19] 

= f^^xt^ , (36) 

where k, and p are dimensionless constants whose values depend on the specific process 
involving the specific kind of TD under consideration we see that Eq. (35) corresponds to 
p = 1 and 

^ ^ ^ (lo^)"' (-"^r (Jk)^ {j^f . (-) 

Using Eqs. (17), (35) and (23), we finally get the following expression for the expected 
fiux in terms of the parameter rj: 

= 3.. 3.^15^)^ (^)^.„.^^^^^^^^^^^ 

_EoY r 1 dE,{Eo,Zi) dN, 2 -2 -1 -1 /gox 

— tt; / ; ; — eV m sec sr , 38 

1020 eVy io l + z, dEo dx ' V ^ 

where x = 2Ei/mx- 

Eigs. 1 and 2 show the expected fiux of various particles for various different injection 
spectra and for different values of nix- The required value of the parameter rj is obtained by 
the normalization procedure described in section 3.5. Eigs. 3 and 4 show the contour plots 
for the normalized value of rj in the Om — parameter space. Erom these plots we see that 
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even if the actual value of Om is significantly less than the Parker bound value, collapse of 
monopolonia can give significant contribution to UHE CR provided rj is sufficiently large (say, 
between 30 and 70). Indeed, as is clear from Eq. (38), the contribution of monopolonia to 
the fiux of cosmic rays increases exponentially with the parameter rj. Of course, for too large 
values of ?y, the monopolonium is so tightly bound that the monopole and the antimonopole 
constituting it do not interact with other monopoles/antimonopoles but are dominated by 
the local mutual coulomb interaction. As a consequence, the equilibrium analysis is not 
expected to be valid for large rj. 

5 Conclusions 

In summary, we have studied the process of monopolonium collapse as a possible source of 
production of the UHE cosmic rays. The properties of the spectra of the UHE nucleons, 
photons and neutrinos produced by this process are essentially same as those already de- 
scribed in earlier works [19, 43] that dealt with UHE cosmic ray production from topological 
defects in general; the reader is advised to consult these references for details of the spectra. 
In the present paper, our main aim has been to estimate the overall magnitude of the fiux 
of UHE cosmic rays produced by monopole annihilations. This is determined not only by 
the actual abundance of monopoles in the universe (which, unfortunately, is unknown), but 
also by the monopolonium abundance relative to monopoles. We have estimated the relative 
monopolonium abundance that would be required in order to produce sufficient fiux of the 
UHE cosmic rays as observed. Eor a monopole abundance saturating the Parker bound, 
say, a fractional monopolonium abundance of ~ 10~® (with nix ~ 10^^ GeV) can give rise 
to measurable UHE CR fiux and can, in principle, explain the recent HECR events [8, 10]. 
The required fractional monopolonium abundance, however, increases with decreasing value 
of the monopole abundance. Within the context of an analysis that assumes monopoles to 
be in thermal equilibrium at the relevant time of monopolonium formation, we find that 
the required relative monopolonium abundances may well be realizable. Out of equilibrium 
processes like transition to a transient superconducting phase [55, 56] giving rise to transient 
magnetic fiux tubes connecting monopole-antimonopole pairs could well enhance monopolo- 
nium abundances beyond the equilibrium values and compatible with the required numbers 
derived above. A full microscopic treatment is, however, required in order to calculate these 
abundances from first principles. Work along this direction is in progress. 
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QCD motivated injection 
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Figure lA: Observable neutrino (short-dashed hne), gamma ray (long-dashed line) and 
proton (solid line) spectrum produced by monopolonium collapse. An X particle mass of 
nix = 10^^ GeV and the QCD-motivated injection spectra discussed in section 3.2 were used. 
The combined proton and gamma ray flux was normalized at fO^^'''eV to the "extragalactic 
flux component" (thin solid line) [7] fitted to the high energy Fly's Eye data (dots with error 
bars) as described in section 3.5. Also shown (dash-dotted line) is an approximate limit on 
the neutrino fiux determined from the non-detection of deeply penetrating particles by the 
Fly's Fye detector [52, 53]. 
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QCD motivated injection 
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Figure IB: Same as Fig. lA but for an X particle mass mx = 10 GeV. 
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Phenomenologically motivated injection 
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Figure 2A: Same as Fig. lA, except that the phenomenologically motivated injection spec- 
tra discussed in section 3.2 were used. Edges in the gamma ray and neutrino spectrum are 
due to simplifying numerical assumptions. 
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Figure 2B: Same as Fig. IB, except that the phenomenologically motivated injection spectra 
discussed in section 3.2 were used. Edges in the gamma ray and neutrino spectrum are due 
to simphfying numerical assumptions. 
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QCD motivated injection 




Log(m /GeV) 
X 



Figure 3: Contours of the relative monopolonium abundance as defined in Eq. (25) 
(solid lines) and the corresponding ?y-parameter within the equilibrium analysis (dashed 
lines) following from the normalization explained in section 3.5. QCD-motivated injection 
as discussed in section 3.2. was assumed. In the nix — ftnh^ plane the ^j- and ?y-contour 
lines correspond to I, O.I, • • • , I0~^ and 70, 60, • • • , 30, respectively, decreasing from the lower 
right to the upper left. The shaded parameter region is excluded within the modified Parker 
bound [37]. 
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